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We recently identified a set of chro- 
mosome domains that are early 
replicating uniquely in pluripotent cells. 
Their switch from early to late replica- 
tion occurs just prior to germ layer com- 
mitment, associated with a stable form of 
gene silencing that is difficult to reverse. 
Here, we discuss results demonstrating 
that these domains are among the least 
sensitive regions in the genome to global 
digestion by either MNase or restriction 
enzymes. This inaccessible chromatin 
state persists whether these regions are 
in their physically distended early repli- 
cating or compact late replicating con- 
figuration, despite dramatic changes in 
3D chromatin folding and long-range 
chromatin interactions, and despite large 
changes in transcriptional activity. This 
contrasts with the strong correlation 
between early replication, accessibility, 
transcriptional activity and open chro- 
matin configuration that is observed 
genome-wide. We put these results in 
context with findings from other studies 
indicating that many structural (DNA 
sequence) and functional (density and 
activity of replication origins) properties 
of developmentally regulated replication 
timing ("switching") domains resemble 
properties of constitutively late replicat- 
ing domains. This suggests that switch- 
ing domains are a type of late replicating 
domain within which both replication 
timing and transcription are subject to 
unique or additional layers of control not 
experienced by the bulk of the genome. 
We predict that understanding the 
unusual structure of these domains will 
reveal a novel principle of chromosome 
folding. 



Introduction 

DNA replication in mammals is regulated 
in a spatio-temporal manner that sug- 
gests a relationship to higher order chro- 
matin structure in the nucleus. During 
early S phase, sites of DNA replication 
are detected throughout the interiorly 
located euchromatic regions of the nucleus 
while regions at the nuclear periphery are 
mainly replicated late during S phase. 1 
Each genomic segment is replicated at a 
distinct time during S phase, which can 
be differentially regulated in the context 
of development and disease. 2,3 Evidence 
for megabase-sized chromosomal domains 
of coordinate replication timing ("replica- 
tion domains"; RDs) comes initially from 
cytological analyses, 4 though it has been 
difficult to gain further insight due to the 
lack of robust means to analyze this level of 
chromosome structure. For instance, elu- 
cidating what defines boundaries separat- 
ing spatially segregated early and late RDs 
is central to our understanding of nuclear 
structure and function, but it has been 
almost impossible to approach this ques- 
tion without knowing where RD bound- 
aries are localized with respect to DNA 
sequences. However, the situation changed 
dramatically with the recent development 
of genome-wide replication timing assays, 
which have comprehensively mapped RDs 
in various cell types. 2,5,6 Our group has 
profiled replication timing of more than 
40 cell types from mouse and human. 
These data sets are available through 
Replicationdomain.org 7 and a select 
group are available through the UCSC 
ENCODE genome browser. Using this 
method, we and others discovered that 
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Figure 1. A model for chromatin reorganization associated with a subnuclear compartment 
switch. 12 In pluripotent stem cells, domains subjected to EtoL regulation partially unfold and 
move into the early-replicating subnuclear chromatin interaction compartment, where they are 
suppressed from interacting with late-replicating chromatin, even neighboring late-replicating 
domains. During differentiation, these domains retract and acquire a similar spatial compaction 
to surrounding domains. In this configuration, they interact more frequently with late-replicating 
chromatin including neighboring domains, while, at the same time, retaining the identity of their 
self-interacting unit boundaries. FISH images show chromatin organization of the Dppa2/4 do- 
main in ESCs and NPCs. Several FISH probes (red for the neighboring late domains and green for 
the entire early domain) were hybridized and detected simultaneously in ESCs and NPCs nuclei to 
visualize chromatin organization of the Dppa2/4 domain. 



RD structures are reorganized genome- 
wide during embryonic stem cell (ESC) 
differentiation into various cell types, 
with approximately 50% of the genome 
undergoing changes in replication timing 
during development. 2,8 Cell-type specific 
reorganization of RDs is coordinated with 
transcriptional changes and is conserved 
between mouse and human. 2,9 

In the course of these studies, we iden- 
tified a set of domains that switch from 
early to late replicating just prior to germ 
layer commitment, approximately coinci- 
dent with X-inactivation, and are late rep- 
licating in all other analyzed cell types. 8 
We showed that these domains undergo 
an unusually stable form of gene silencing, 
in that both late replication and expression 
of genes within these domains is particu- 
larly resistant to reprogramming back to 
the pluripotent state. 8 Here, we highlight 



a recent study of the large-scale structure 
of these domains that has revealed the 
surprising finding that these domains 
undergo dramatic changes in their 3D 
organization and subnuclear chromatin 
interaction compartments while remain- 
ing exceptionally inaccessible to nucleases. 
We discuss these findings in the context 
of the general properties of developmen- 
tally regulated vs. constitutively replicated 
domains, and the extent to which this may 
be a more widespread mechanism. 

A Novel Principle 
of Chromatin Folding 

Domains that replicate early uniquely in 
the pluripotent cell state represent approx- 
imately 6% of the genome. 8,10 What is 
the biological significance of these early 
to late (EtoL) domains? Both reversal of 



directed differentiation experiments and 
a comparison of genome-wide replication 
timing between fully reprogrammed and 
partially reprogrammed induced pluripo- 
tent stem cells (iPSCs) identified these 
EtoL switching domains as regions that 
resist replication-timing as well as tran- 
scriptional reprogramming, coincident 
with their switch to late replication. 8 Since 
these EtoL domains contain pluripotency 
associated genes such as Rexl and Dppa2/4 
that occlude the binding of pluripotent 
cell transcription factors during the early 
stages of the reprogramming process, 11 
late replication may be linked to the stable 
silencing of these pluripotency-associated 
genes. What mechanism could constitute 
this link? One possibility is that EtoL 
switching domains form a specific chro- 
matin structure that is particularly diffi- 
cult to reprogram. To test this hypothesis, 
we examined chromatin accessibility in 
replication timing switching and non- 
switching domains genome-wide during 
differentiation of ESCs to neural precursor 
cells (NPCs) and found that EtoL switch- 
ing domains do have unusual chromatin 
accessibility pattern. 12 First, MNase was 
used to digest chromatin from fixed nuclei 
in order to make direct comparisons to the 
structure of these domains determined by 
fluorescence in situ hybridization (FISH). 
We confirmed a strong genome-wide cor- 
relation of general accessibility to early 
replication previously shown by others, 13 ' 14 
but found that EtoL switching domains, 
including Dppa2/4 and Rexl domains, 
formed inaccessible chromatin before and 
after the replication-timing switch. In 
fact, we could detect almost no changes in 
accessibility despite large swings in replica- 
tion timing. The same result was obtained 
when unfixed chromatin was used for 
MNase digestion. To confirm these results 
with higher molecular weight nucleases, 
we also performed digestions using a com- 
bination of restriction enzymes (Nlalll 
and Mspl). To our knowledge, this is 
the first time restriction enzymes have 
been employed to probe global chroma- 
tin accessibility, but we obtained identical 
results to MNase digestion. 

The fact that early replicating Dppa2/4 
and Rexl domains in ESCs form inacces- 
sible chromatin but harbor transcription- 
ally active genes within them raises the 
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possibility that another layer of chromatin 
structural change is occurring to form the 
transcriptionally permissive chromatin 
domain. Indeed, when we examined the 
distance between pairs of FISH probes 
spanning across these domains in the 
nucleus, we found an ESC-specific dra- 
matic chromatin unfolding restricted to 
the region of replication timing switch, 
which returned to the normal degree of 
chromatin condensation after differentia- 
tion (Fig. 1). Observed chromatin reorga- 
nization at the cytological level suggests a 
spatial eviction of early Dppa2l4 and Rexl 
domains from neighboring late replicating 
domains in ESCs and led us to examine 
intranuclear chromatin interaction pat- 
terns of these domains before and after 
cell differentiation. We previously dem- 
onstrated that replication timing strongly 
correlates with alternate compartments of 
chromatin interactions defined by Hi-C 
(high resolution chromatin conforma- 
tion capture). 9 ' 15 Hence, we performed 4C 
(chromosome conformation capture on 
chip) using sites centrally located within 
the Dppa2l4 and Rexl domains as baits. 
We found that, when early replicating, 
these domains interact with early repli- 
cating chromatin, while after switching 
to late replication, they interact with late 
replicating chromatin, providing the first 
demonstration of a dynamic switch in 
interaction compartments induced by dif- 
ferentiation. These findings were strength- 
ened by a recent Hi-C study showing 
spatial partitioning of Dppa2l4 and Rexl 
domains from neighboring chromosomal 
domains in mouse ESCs but not in neu- 
ral cortex. 16 Taken together, the data show 
that a chromatin interaction compart- 
ment switch — without a change in acces- 
sibility — distinguishes these chromatin 
domains before and after differentiation, 
unveiling a novel principle of chromatin 
folding at developmentally regulated RDs. 

Unfolded Inaccessible Chromatin: 
No Longer an Oxymoron 

Unfolding chromatin and moving it to a 
different sub-nuclear chromatin interac- 
tion compartment may also explain how 
genes in closed chromatin domains can be 
transcribed. 13 But why do two antagoniz- 
ing mechanisms for chromatin regulation 



(unfolded yet inaccessible chromatin) co- 
exist in chromosomal domains containing 
pluripotency-associated genes? Given that 
these genes are expressed almost exclu- 
sively in pluripotent stem cells but remain 
silenced throughout the rest of develop- 
ment, it is possible that genes in closed 
chromatin are readily silenced in response 
to differentiation signals and are prevented 
from aberrant activation. Previous studies 
have revealed similar large-scale chroma- 
tin unfolding events, but some of these 
are closely linked to transcription, 17 " 19 
while others appear to be independent 
of transcription. 20,21 This suggests that 
there are multiple mechanisms leading to 
large scale unfolding, only some of which 
require active transcription. Our RNA- 
DNA FISH and transcription inhibition 
experiments demonstrated that active 
transcription is not necessary to maintain 
the ESC-specific chromatin configura- 
tion at Dppa2l4 or Rexl domains, 12 sug- 
gesting that the distended chromatin state 
described here reflects transcriptional 
competence of these domains rather than 
transcription per se. 

Chromatin Interaction 
Compartments 
and Transcriptional Competence 

The nucleus is highly compartmental- 
ized into sub-nuclear domains in which 
specific sets of proteins are concentrated 
for various nuclear processes. It is increas- 
ingly recognized that interaction between 
these nuclear domains and chromatin, at 
least in part, contributes to establish and 
maintain cell type specific nuclear archi- 
tecture and function. 22,23 We observed 
a chromatin interaction compartment 
switch accompanied by a replication tim- 
ing switch at developmentally regulated 
RDs. We still have no clear answer for 
what this compartment switch means, but 
several lines of evidence suggest a link to 
transcriptional competence. It has been 
shown that actively transcribed genes 
preferentially colocalize at "transcription 
factories" in which RNA polymerases 
are highly concentrated. 24 Transcription 
factories are found in the interior of the 
nucleus and close proximity of these facto- 
ries to early replicating chromatin suggests 
a functional coupling of these two events 



in the nucleus. 25 This close relationship is 
supported by genome-wide studies show- 
ing that early replication is strongly cor- 
related with transcriptional activity. 2,26 " 28 
On the other hand, a lot of late replicat- 
ing chromatin is localized near the nuclear 
periphery where transcriptional activity is 
generally low. 29,30 In some gene loci, local- 
ization to the nuclear periphery is associ- 
ated with transcriptional repression. 30 " 34 
Tethering of endogenous gene loci to the 
nuclear periphery can induce downregula- 
tion of some active genes. 30,35,36 Given the 
close relationship between late replication 
and silent nuclear domains, preferential 
interaction with neighboring late replicat- 
ing chromatin compartment observed at 
Dppa2/4 and Rexl domains may increase 
interaction of these chromosomal domains 
with transcriptionally repressive nuclear 
compartments. Indeed, Dppa2/4 and Rexl 
form lamina-associated domains (LADs), 
transcriptionally repressive perinuclear 
domains, upon ESC differentiation to 
NPCs. 37 Interestingly, although LADs 
are smaller and do not change as much 
during differentiation as replication tim- 
ing, LAD and RD boundaries align very 
closely at the Dppa2/4 and Rexl domains 
(Fig. 2). What mechanism could link sub- 
nuclear compartments, replication timing 
and transcriptional competence? It was 
reported that extra-chromosomally repli- 
cating viral DNA can switch from acety- 
lated to de-acetylated chromatin when 
stochastically replicating either early or 
late during consecutive S phases. 38 Hence 
one attractive hypothesis is that early and 
late replicating compartments are enriched 
with different sets of chromatin modifiers 
and therefore facilitate the formation of 
transcriptionally permissive and repres- 
sive chromatin structure, respectively. To 
date, there is little evidence for S-phase 
specific subnuclear compartmentaliza- 
tion of chromatin modifiers, 1 although 
that HDAC3 has shown to interact with 
Lap2b and localize specifically to nuclear 
periphery. 30,39 

New Dimensions 
to DNA Replication Research 

What is regulating RD structure? The 
conservation of replication timing and RD 
structures between mouse and human. 9,40 
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Figure 2. Replication timing (A) and Lamin B1 binding (B) of the Dppa2/4 domain. The replication 
timing-switching domain is indicated in the red box. Figure was assembled using the Replication- 
Domain.org database. 7 



and the consequences of chromosome 
rearrangements within RDs, 3,41 suggest 
the existence of cis-regulatory elements 
that define these structural properties. On 
the other hand, dynamic reorganization of 
RD structure during cell differentiation 
suggests that epigenetic mechanisms also 
take part in the formation of RD structure. 
We previously showed that H3K9me2 is 
highly enriched at late replicating periph- 
eral chromatin. 42 Surprisingly, RD struc- 
tures are not affected in cells lacking G9a, 
the histone H3K9 dimethylase, despite 
the fact that mutant cells display loss of 
H3K9me2 mark and aberrant expression 
of genes within late replicating peripheral 
chromatin. 43 It has also been reported that 
mutations in other chromatin and DNA 
modifiers can show profound effects on 



gene expression patterns genome-wide 
but have little or no effect on replication 
timing. 44 By contrast, the only gene dis- 
ruption shown to affect the global repli- 
cation-timing program is an enigmatic 
nuclear matrix protein called Rifl, disrup- 
tion of which has subtle effects on tran- 
scription. 45,46 These results suggest that 
neither the structural changes in chroma- 
tin caused by these mutations nor the sub- 
sequent changes in transcriptional activity 
are sufficient to drive this level of chromo- 
some reorganization. 

RD boundaries closely align to topo- 
logical chromosomal domain boundar- 
ies, 912 ' 16 implying that the boundaries 
between units of chromatin that can self- 
interact also form the boundaries of units 
of replication regulation. Considering that 



chromatin looping by long-range interac- 
tion is a major mechanism behind topo- 
logical domain formation, it is possible 
that factors mediating chromatin loop for- 
mation may play some roles in the regula- 
tion of RD structure. CTCF is one of the 
well-known proteins mediating chromatin 
loop formation. It has been shown that 
CTCF is enriched in topological domain 
boundaries revealed by Hi-C. 16 The obser- 
vation that only 15% among all CTCF 
binding sites are associated with these 
boundaries suggests either that CTCF 
may require additional factors to confer 
its specific action at boundaries or that 
CTCF binding may reflect occupancy of 
some other regulatory factors at boundar- 
ies. Currently, little is known about roles 
of CTCF and its interacting proteins in 
RD formation. Interestingly, knockdown 
of Cohesin induces inefficient origin acti- 
vation accompanied by increased chro- 
matin loop size. 47 Although Cohesin and 
CTCF are known to cooperatively main- 
tain chromatin loops, observed effects of 
Cohesin deficiency are independent of 
CTCF. Recently, it was reported that rep- 
lication timing is affected by the binding 
of special AT-rich sequence binding pro- 
tein 1 (SATB1) which is know to mediate 
chromatin loop formation. 48 Importantly 
however, while many proteins have shown 
to be involved in chromatin loop forma- 
tion in mammalian cells, their roles in 
RD formation have not been examined 
and many loops that are formed and 
broken during transcriptional regulation 
have no effect on chromatin interaction 
compartments . 45 

A study in yeast provides a detailed 
insight into the link between RD orga- 
nization and long-range chromatin inter- 
action. Yeast transcription factors Fkhl 
and Fkh2 have been shown to regulate 
global replication timing potentially by 
bringing origins in close proximity in the 
nucleus. 50 These proteins bind near origins 
and also interact with origin recognition 
complex (ORC), which possibly allows 
them to bridge between origins to form a 
cluster. Replication timing is determined 
by the association of origins with limit- 
ing amounts of replication initiation fac- 
tors such as Cdc45 during G [ phase. 51 " 53 
As loss of Fkhl and Fkh2 results in 
either advanced or delayed origin firing, 
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it is possible that origin clustering could 
affect spatial proximity of origins to rep- 
lication initiation factors-enriched nuclear 
domains, thereby regulating their associa- 
tion during Gj phase. Recently Rifl has 
been shown to be required for the global 
replication-timing program both in yeast 
and mammals. 45,46,54 In Rifl -deficient cells, 
some regions undergo EtoL replication 
timing switching, while others undergo 
LtoE switching. Perturbed spatial dis- 
tribution of replication foci patterns and 
increased chromatin loop size are observed 
in 7?zy?-deficient mammalian cells. It will 
be interesting to examine effect of Rifl 
mutation on chromatin interaction land- 
scape in the nucleus, which may provide 
additional insight into the link between 
RD organization and long-range chroma- 
tin interaction. Thus, a topological view 
of DNA replication control opens new 
avenues for DNA replication research. 

Switching Domains are a Type 
of Late Replicating Domain 

Our finding that developmentally regu- 
lated RDs (switching domains) remain 
inaccessible even when early replicat- 
ing supports our earlier contention that 
switching domains share properties of late 
replicating domains, as if late replication is 
the default replication timing for switch- 
ing domains. 55 Indeed, many observations 
now suggest that there are general rules to 
set the default replication program, but 
that some additional property must be 
imposed upon switching domains for them 
to replicate early. First, the DNA sequence 
composition of switching regions, includ- 
ing GC content and the density of LINE 
elements, are similar between late and 
switching domains; for example, more 
than 6% of the sequence composition of 
both late and switching regions consist 
of LINE elements (Fig. 3). 55 Functional 
properties other than chromatin acces- 
sibility are also shared between late and 
switching regions. A recent genome-wide 
origin mapping study has shown that, as 
expected, early replication is positively 
correlated with the density and the activ- 
ity (frequency of usage) of replication ori- 
gins but, unexpectedly, switching regions 
have low origin density and activity even 
when they are early replicating. 56,57 This 
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Figure 3. Replication timing changes predominantly in LINE-rich regions. Replication timing 
values (log 2 early/late-S enrichment) are plotted for RefSeq genes in ESCs and NPCs, with values 
in the lower right and upper left representing LtoE and EtoL timing changes, respectively. Points 
are colored by the percentage of LINE-1 sequence assigned to each gene from the nearest 100 kb 
window, from blue (> 6% LINE) to yellow (< 6%LINE). 
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Figure 4. Timing and accessibility are uncoupled in replication timing switching regions. Correla- 
tion between replication timing values and accessibility levels (log 2 nuclease accessible/inacces- 
sible fractions) in ESCs and NPCs, ranked by the level of developmental timing changes across 9 
mouse cell types. Correlations were measured in an overlapping moving window of 400/10,566 
200 kb regions ranked by most static to most dynamic replication timing (smallest to largest 
standard deviation in timing values) during development. 
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observation is consistent with a model in 
which general regulation of replication 
timing within switching domains renders 
them late replicating, while early replica- 
tion is mediated by some mechanism that 
can activate the relatively low number of 
potential origins within these domains 
earlier during S phase, without increasing 
the efficiency of their usage. Alternatively, 
a generalized increase in origin usage 
could remain undetected by the methods 
used in this report if there was an increase 
in low efficiency origin activity dispersed 
throughout a much larger portion of 
the domains. Regardless of the mecha- 
nism, these observations together suggest 
that developmentally regulated RDs are 
an independent breed of chromosomal 
domain with unusual regulatory proper- 
ties. All the while, there are some proper- 
ties of chromatin that do globally correlate 
with replication timing in the switching 
domains to an equal extent as the bulk of 
the genome, such as correlations to tran- 
scription and chromatin post-translational 
modifications associated with transcrip- 
tional activity.' Understanding what sep- 
arates these properties that distinguish 
switching domains from constitutively 
late replicating domains should provide 
clues to the nature of this breed of chro- 
mosomal domain. 

How Widespread Is This Novel 
Principle of Chromatin Folding? 

The fact that 50% of the genome 
switches replication timing at some stage 
of development (possibly more will be 
revealed when other tissues, including 
extra-embryonic tissues, are explored) 
implies that if this novel mechanism 
applies to all switching regions, it is not 
unique but quite widespread. At pres- 
ent, the answer to this question is not 
clear. Based on replication timing pro- 
files of 9 mouse cell types (Fig. 4), 
it is clear that accessibility and replica- 
tion timing become progressively uncou- 
pled the more dynamic the replication 
time of a domain is during development. 
However, FISH studies of 3D folding 
are laborious and single-locus based and 
it is not yet clear whether either Hi-C 
or replication timing is a reliable proxy 
for the scale of 3D folding that we have 
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Figure 5. Low accessibility is a shared property of replication timing switching genes. (A) Acces- 
sibility level (log 2 nuclease accessible/inaccessible fractions) as a function of developmental rep- 
lication timing changes between ESCs and NPCs. In contrast to the general relationship between 
accessibility and early replication, timing switching regions in both directions (EtoL or LtoE) are 
relatively inaccessible even in ESCs (particularly EtoL), prior to observed changes in subnuclear 
organization and chromosomal contacts at EtoL loci. (B) Accessibility vs. median absolute devel- 
opmental timing changes across all cell types. Regions with large or frequent timing changes in 
mouse development predict low accessibility in ESCs and NPCs, despite these changes arising in 
different lineages. For both graphs, accessibility and timing values are measured at RefSeq gene 
promoters, with a shaded 95% confidence interval about the mean. 
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reported. Moreover, Hi-C maps are still 
very expensive to generate. While maps 
of DNasel hypersensitive sites are now 
available for many mouse and human cell 
types, 58,59 maps of global chromatin acces- 
sibility are limited. Also, our study indi- 
cated that while switching domains were 
generally less accessible than the bulk of 
the genome, regions switching LtoE dur- 
ing ESC differentiation were considerably 
more accessible than those switching EtoL 
(Fig. 5A). Since accessibility data were 
only collected during one single cell fate 
change, it is difficult to predict how RDs 
switching "EtoL" or "LtoE" behave dur- 
ing other cell fate transitions. Moreover, 
not all replication timing switch domains 
resist reprogramming when differen- 
tiation is reversed. 8 Nonetheless, regions 
with large or frequent timing changes 
during mouse development predict low 
accessibility in ESCs and NPCs, despite 
these changes arising in different lin- 
eages (Fig. 5B). Hence, while we certainly 
expect many more domains in the genome 
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